Momentum Distribution and Renormalization Factor in Sodium and the Electron Gas 
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We present experimental and theoretical results on the momentum distribution and the quasipar- 
ticle renormalization factor in sodium. From an x-ray Compton-profile measurement of the valence- 
electron momentum density, we derive its discontinuity at the Fermi wavevector. This yields an 
accurate measure of the renormalization factor that we compare with quantum Monte Carlo and 
Go Wo calculations performed both on crystalline sodium and on the homogeneous electron gas. Our 
calculated results are in good agreement with the experiment. 

PACS numbers: 71.10.-w, 78.70.Ck, 71.20.Dg, 02.70.Ss 
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Introduction The homogeneous electron gas (HEG), 
also known as jcllium, is one of the most fundamental 
models in condensed matter physics It is one of the 
simplest many-body systems which can still describe sev- 
eral properties of real solids, especially of the alkali met- 
als. For almost half a century, the accurate description 
of many-body correlation effects has challenged quantum 
many-body theory and HEG is the canonical workbench 
to test different theoretical methods 0-11]. Although 
the exact analytic solution of the many-body problem 
in HEG is still unknown, today quantum Monte Carlo 
(QMC) calculations are widely accepted to provide the 
most reliable results on e.g. the correlation energy. The 
situation is less clear concerning spectroscopic quantities 
such as the momentum distribution, n(p). The accu- 
racy of the theoretical methods in this respect is not well 
understood, different approaches yielding a wide range of 
varying results. This fundamental issue has remained un- 
resolved, mainly due to a lack of accurate, bulk-sensitive 
and unambiguous experimental probes that could be used 
to compare the theories with. 

Experimentally, one of nature's closest realizations of 
HEG is formed by the valence electrons in alkali metals, 
especially Na. Here, we present very accurate experi- 
mental and theoretical results on the electron-momentum 
distribution of Na. The single occupied valence band of 
Na has an almost spherical Fermi surface and its prop- 
erties in ambient conditions with a density parameter 
r s =3.99 can be directly compared with theoretical re- 
sults on HEG. In particular, we obtain a precise experi- 
mental reference value for the quasiparticle renormaliza- 
tion factor, Zk F , which characterizes the discontinuity of 
the momentum distribution at the Fermi surface at this 
density 

From the Compton profile (CP) measured by inelastic 
x-ray scattering experiments on bulk sodium, we derive 
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FIG. 1: (Color online) The momentum distribution of Na 
determined by experiment, QMC using Slater- Jastrow wave- 
functions, GoWo, and DFT-LDA calculations. The ideal- 
Fermi gas step function is also shown. 



n(p) and obtain = 0.58(7). Compared to previous 
experiments our experimental resolution provides a 
clear observation of the discontinuity at the Fermi surface 
in a direct and mo del- independent way. 

We compare our experimental results to theoreti- 
cal calculations using QMC and Go Wo methods, both 
done for HEG and for solid Na taking into account 
the electron-electron interaction and band-structure ef- 
fects. Our calculations confirm the jellium-likc behavior 
of Na and allow us to quantify the small band-structure- 
induced deviations from HEG. Finally, we compare the 
results with other many-body approximations applied to 
HEG in literature [|H1] (Table Unless explicitly spec- 
ified, we use atomic units (a.u.). 

Theory The momentum distribution per spin state, 
n(p) (see Fig. 1), is one of the basic many-body observ- 
ables where the Pauli principle for fermions is directly 
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0.57(7) 0.58(7) 
0.68(2) 0.70(2) 0.69(1) 
0.66(2) 

0.64(1) 0.65(1) 0.64 [2] 
0.793 [6] 
0.45 [4, 5j 
0.59 [4]~ 
0.61 [8J 
0.615 [5J 

0.71 [7] 



TABLE I: Our results for ( and Z kp of Na and HEG at 
r s = 4.0: experimental value compared to GoWo, as well 
as QMC results based on Slater-Jastrow (QMC-SJ) and more 
precise backflow (QMC-BF) wavefunctions. Various theoret- 
ical values of HEG from literature are also quoted. 



visible. It is the probability to observe an electron with 
momentum p. For a non-interacting Fermi gas at zero 
temperature, n{p) is 1 for p below the Fermi momentum 
Pf and above, i.e. n(p) = 9(pf —p) with a discontinuity 
C = n{Pp) — n(p~t) = 1 occurring at the Fermi surface. 
For a non-interacting crystalline system, the electrons 
occupy Bloch wavefunctions, 0„k(r) = J2g ^Sc e ^ k+G ^ r i 
where k is the crystal momentum, v the band index, and 
G are reciprocal lattice vectors. For systems like Na with 
one valence band (y = 1) whose Fermi surface is entirely 
contained within the first Brillouin zone (1BZ), the band 
structure reduces the discontinuity of the momentum dis- 
tribution, £ = l0^Ti°k | 2 < 1- From a density- functional 
theory calculation within the local-density approxima- 
tion (DFT-LDA) for Na, we obtain Cdf T = 0.98(1). The 
calculated valence band is an almost perfect parabola, 
its wavcfunction is nearly isotropic, and its Fermi sur- 
face deviates from a perfect sphere by only 0.2%. These 
deviations of the Fermi surface from a perfect sphere nec- 
essarily lead to a further, albeit small, reduction of the 
discontinuity when n(p) is orientationally averaged. 

Many-body effects introduce a much larger reduction 
of the discontinuity at the Fermi surface. This is known 
as the quasiparticle rcnormalization factor, Z^ F . The- 
oretical predictions for Z^ F using different approxima- 
tions range from 0.45 to 0.79 for the density consid- 
ered (see Table I). In general, the renormalization factor 
Zk F is related to the self-energy E^^fcjw) via Z^ F = 
(1 - <9Ei,i (k^ 1 , w)/ dui\ UJ =£ F ) — 1 , E = E c " c + S C "P con- 
taming the electron-electron interactions, S c_c , and the 
electron-phonon effects, S c_p . For HEG the discontinu- 
ity in the momentum distribution is £heg = Z^ F . In 
a jellium-like system such as Na, band-structure effects 
and many-body correlations can be factorized so that 
CNa = l0^Ti°k F | 2 ^k F with a renormalization factor very 
close to the value for HEG at the same density, if phonon 
effects can be neglected. 

To determine Z^ F theoretically, we performed pscu- 
dopotential diffusion QMC 1(| ll| calculations of bulk 



sodium based on a Slater-Jastrow (SJ) wavcfunction us- 
ing the QMCPACK code, and more precise calculations 
using backflow (BF) for HEG. Complementary to QMC, 
we have done a non-self-consistent (one-shot) Go Wo cal- 
culation Q starting from the DFT-LDA electronic struc- 
ture using the ABINIT code. 

Within both methods, pseudopotentials are used to 
describe the core electrons, based on a regular static 
lattice for the ions, neglecting effects due to electron- 
phonon coupling. Whereas core correlation effects only 
give smooth corrections that do not influence the value 
of the renormalization factor, electron-phonon coupling 
may lead to a further decrease of Z^ F . However, since 
the phonon Debye frequency uid is small compared to 
the Fermi energy, main effects of £ e " p are expected only 
within a narrow momentum region around pp, with 
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As we will see below, those 



effects are beyond the resolution of the experiment. 
The static approximation and the use of pseudopotcn- 
tial should thus be sufficient to obtain the value of Z-^ F , 
whereas they may be less accurate to predict the whole 
CP. 

Experiment A unique bulk-sensitive probe of the mo- 
mentum density is offered by Compton scattering of x- 
rays [l2T ] . The experiment measures the spectra of x-rays 
scattered by an electron system. When the energy trans- 
ferred to the electron is much larger than its binding 
energy, the so-called impulse approximation (IA) is valid 
and the measured spectrum is related to the CP, which 
in isotropic average normalized to one electron is 



8tt Pf 



dtt p n(p) dp. 



(1) 



Here q is the component of the ground-state momentum 
of the electron projected onto the scattering vector. As- 
suming an isotropic system, n(p) can thus be extracted 
by a differentiation of the CP, 



n{p) 



2p% dJ(q) 



'dp dq 



q=p 



(2) 



For the non-interacting HEG the CP is an inverted 
parabola J(q) = -^pr(pp — q 2 ) for q < pf and vanishes 
for q > Pf {pf = 0.49(1) a.u. for Na). Many-body ef- 
fects promote a part of the electrons from below to above 
Pf- The CP, while being always continuous, should re- 
tain a kink, i.e. a discontinuity in the first derivative at 
ftF- Measuring the CP accurately allows the extraction 
of n(p) by using Eq. ([2]). The determination of Zk F via x- 
ray Compton scattering has been a long-standing goal of 
many scientists jl2j ]. The simultaneous requirements of 
extremely high momentum resolution and statistical ac- 
curacy as well as difficulties in separating band-structure 
and correlation effects have made such attempts difficult, 
leading to anomalously low values of e.g. Z^ 1 = 0.1(1) 
for Li (r s = 3.25) [13]. Promising results were given 
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FIG. 2: (Color online) The measured x-ray-scattering spectra 
from Na as a function of energy transfer, for both experimen- 
tal runs. The experimental spectra consist of overlapping va- 
lence and core contributions. Theoretical core contributions 
are shown for both QSCF and FEFFq treatments. 



for Al (r s = 2.07) [14| by comparisons with an analyt- 
ical model of n(p) with an adjustable Z^ F giving 
the best agreement with Z^ 1 m 0.7. This determination 
however assumed a specific shape of n(p) and thus was 
not model- independent. Our choice of a HEG-like system 
of Na combined with ultra-high resolution measurements 



allows us to accurately determine the n(p) and Z^ in a 
modcl-indcpcndcnt way. 

The experiments were performed at the bcamlinc ID16 
fl5j of the European Synchrotron Radiation Facility on 
polycrystalline sodium. The spectrometer was based on 
a Rowland circle with spherically bent analyzer crystals 
with a Bragg angle of 89°. The measurements were done 
by changing the incident-photon energy E\ and observ- 
ing the flux of scattered photons into a fixed scattering 
angle 29 at a fixed energy E 2 ■ We used two different con- 
figurations and photon-energy ranges to verify the result 
in two independent ways. In the first experiment (run 1), 
we used a single Si(555) analyzer crystal, E 2 = 9.9 keV, 
29 = 147° and E 1 = 9.9-11.0 keV. In the second ex- 
periment (run 2), two Si(880) analyzers were used, with 
E 2 = 12.9 keV, 29 = 149° and E x = 12.9-14.0 keV. The 
sample was prepared in a glove box and transported to 
the beamline within an argon atmosphere and pumped 
into a vacuum of 10 -6 mbar. There was no observable 
degradation of the sample when it was inspected after 
the experiment. The measured signal was corrected for 
sample self-absorption as well as changes in the incident 
photon flux and the spectra were measured repeatedly to 
identify any possible instabilities during the experiment. 
None were found and the spectra were finally averaged. 

The measured spectra as a function of energy transfer 
are shown in Fig. [2] The Na L edges are seen at 30-60 
eV, and K edge at 1.07 keV. In run 1 the CP is cen- 
tered at an energy transfer of 365 eV and in run 2 at 
645 eV. Since our interest is in the valence-electron CP, 



FIG. 3: (Color online) The experimentally determined valence 
CPs from the two runs (averaged over q < and q > 0), com- 
pared to the results of Eisenberger et al. [i| , and to the QMC 
CP in the (100) direction. Inset: zoom on the discontinuity 
(points) and the linear fits to the CP around the Fermi mo- 
mentum. Both runs give the same result within the errorbars, 
the fits shown here are from run 2. 



the core contribution has to be subtracted first. Since 
the IA is not strictly valid for the core-electron spectra 
in these experiments, we calculated them with two inde- 
pendent methods: using i) the quasi-self-consistcnt field 
(QSCF) approximation |l6j and ii) the real-space multi- 
ple scattering approach with the FEFF program 17j with 
modifications for calculating the momentum-dependent 

The differences between the 
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scattering cross-section 
two approaches are negligible. The core contribution can 
then be reliably subtracted from the experimental spec- 
tra. The spectra can now be converted into the CP 19| : 
for each energy transfer we can evaluate the scattering- 
electron momentum component q and the measured in- 
tensity is related to the probability of finding the electron 
with that q. 

A finite experimental accuracy in the determination 
of q will introduce a broadening of any sharp features 
in the experimental data. This uncertainty is caused in 
the present experiments by the spread of scattering an- 
gles of the detected radiation. This geometrical contri- 
bution to the g-resolution was Aq = 0.018 a.u. (run 1) 
and Aq = 0.027 a.u. (run 2) full-width-at-half-maximum 
(FWHM). Final-state effects (HE!, i.e. the interaction 
of the scattering electron and the rest of the electron gas, 
are known to cause further broadening of the measured 
valence Compton profiles. We calculated the magnitude 
of this broadening [23], and found it to be effectively 
an additional Gaussian smoothing of 0.08 a.u. (run 1), 
and 0.03 a.u. (run 2) (FWHM). This combined with the 
geometrical resolutions yields effective experimental q- 
resolutions of 0.08 a.u. (run 1) and 0.04 a.u. (run 2). 

Results and discussion The result of the experiment, 
after the analysis described above, is the valence CP 
shown in Fig. [3] The valence CP of a real metal in general 
deviates from the jcllium parabola due to two reasons: (i) 
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correlation modifies the n(p) introducing tails for p > pp, 
and (ii) electron- ion interaction modifies the overall wave- 
function and induces tails for p > pp due to core- 
orthogonalization and the high-momentum components 
<f>rff°- As discussed above, the valence electron wave- 
function of Na is fully contained inside 1BZ and is highly 
free-electron-like, with negligible high-momentum com- 
ponents. The band structure only leads to small (< 3%) 
lowering of the momentum distribution for p < pp, as 
can be seen in the difference between the ideal-Fermi-gas 
and the DFT-LDA results in Fig. 1 . For this reason, we 
can compare the experimental momentum distribution 
and the CP to those of HEG after taking these small 
corrections into account. 

In Fig. [31 the Fermi momentum can be directly seen 
as the discontinuity of the valence CP derivative. From 
the experimental data we deduce pp = 0.49(1) a.u. (LDA 
value 0.481 a.u.). The best determination of the magni- 
tude of the discontinuity at the Fermi surface is provided 
by linear fits to the measured points in the immediate 
vicinity of pp. An inset to Fig. [3] shows these fits, here for 
both negative and positive q from the higher-g-resolution 
run 2. The difference of the slopes for the two sides gives 
C~ = 0.59(7) and ( + = 0.55(7), which allows us to quote 
the average value as = 0.57(7). The errorbar is 
based on the statistical noise of J(q) and the uncertainty 
of pp (cf. Eq. (2)). Using the pure band-structure value 
in LDA, |^J~i k F | 2 = 0.98, we deduce the experimental 
Zk F = Cc N xp/lte°k F | 2 = 0.58(7). We simulated the ef- 
fect of the finite experimental g-resolution by convoluting 
the QMC CP with the resolution function of run 2, and 
repeating the analysis described above. The result was 
an effective lowering of Cqmc 0.02. Since this effect is 
smaller than our errorbar of 0.07, no further correction 
to the final experimental result was found necessary. 

The momentum distribution can be calculated from 
the experimental CP using Eq. (2) and is shown in Fig. 
1. In the differentiation, the effect of statistical noise 
increases and thus we only show the result after averaging 
the values of p < and p > 0, as well as averaging 
adjacent measured points. Due to this, the quantitative 
determination of £ is better done by directly analyzing 
the CP as described above. However, the trend of Table 
I is clearly seen also in Fig. 1 and the experimental n(j>) 
is consistent with the obtained C^p- 

We have also calculated the CP by QMC in the (100) 
direction (Fig. [3]); small differences with respect to the 
experiment remain. Further calculations are necessary to 
study if the inclusion of core electrons or phonon effects 
are necessary to reduce the residual discrepancy between 
theory and experiment in the CP. In Fig. [TJ we show the 
resulting direction-averaged n(p). Finite-size corrections 
[22j are important to determine the structure around pp, 
in particular to obtain the jump at the Fermi surface, 
C£ c = 0.68(2). Within the QMC approach, we can 



determine bare band-structure effects by turning off the 
explicit electron-electron correlations in the underlying 
many-body wavefunction which yields a band-structure 
contribution |<^Ti°k | 2 = 0.97(1) compatible with that 
of DFT-LDA. We therefore obtain = 0.70(2), in 

agreement with QMC calculations of HEG using the same 
type of wavefunctions (SJ). Within HEG, we have per- 
formed calculations using more accurate Slater- Jastrow 
backflow (BF) wavefunctions which indicate a slightly 
lower value, Zff G = 0.66(2). 

From the Go Wo self-energy, we directly obtain Z^ F = 
(1 - d'Ei A (h F ,oj)/duj\ u;=eF )- 1 = 0.65(1). In order to 
determine the jump in n(p), we further need the quasi- 
particle weight which, within GqWq, coincides with the 
Kohn-Sham orbital of the DFT-LDA calculation at the 
Fermi energy, and we get CgoW = 0-64(1). The momen- 
tum distribution within Go Wo is very close to the QMC 
result. 

The agreement between QMC-BF and G0W0 is re- 
markable. Furthermore, within both theories band- 
structure effects can be factorized from correlations 
within the accuracy of the calculation, which enables a di- 
rect comparison with HEG. In Table I, we summarize our 
experimental and theoretical results on Z\^ F and compare 
them to various theoretical results obtained for HEG. On 
theoretical grounds, QMC-BF is considered to give the 
most precise result. Together with Go Wo, it is in reason- 
able agreement with the experiment. 

Our experimental and theoretical values clearly ex- 
clude two different classes of approximation, and thus 
resolve a long-standing theoretical controversy. Whereas 
the so-called on-shell approximation of the RPA 0, S| 
leads to an underestimation of Z^ F , fully self-consistent 
GW calculations Q overestimate its value. The latter 
result shows in particular that the use of the theoret- 
ically more appealing conserving approximation (GW) 
does not result in an improved description of spectral 
quantities compared to non-self-consistent (Go Wo) treat- 
ments. 

Conclusions We have determined the momentum dis- 
tribution of Na valence electrons experimentally and the- 
oretically. In particular, we have related the discontinu- 
ity at the Fermi level to the quasiparticle rcnormalization 
factor of HEG at r s — 3.99, giving a long-sought refer- 
ence value for this fundamental quantity. 
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